# calculate voigt function
from scipy import integrate

def voigtslow(a,u):
    # in tegral of this gives the voigt function
    return exp(-y*y)/((u-y)**2 +a*a)

y = np.linspace(-5,5,100)
a = 0.01
u = 1
plot(y*a+u,f(a,u,y*a+u),'-')

# want to integrate this from -oo to +oo
